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Abstract 

The appearance of quark phases in the dense interior of neutron stars 
provides one possibility to soften the equation of state (EOS) of neutron 
star matter at high densities. This softening leads to more compact 
equilibrium configurations of neutron stars compared to pure hadronic 
stars of the same mass. We investigate the question to which amount 
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^Heisenberg Fellow. 
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the compactness of a neutron star can be attributed to the presence of 
a quark phase. For this purpose we employ several hadronic EOS in the 
framework of the relativistic mean-field (RMF) model and an extended 
MIT bag model to describe the quark phase. We find that - almost 
independent of the model parameters - the radius of a pure hadronic 
neutron star gets typically reduced by 20 — 30% if a pure quark phase in 
the center of the star does exist. For some EOS we furthermore find the 
possibility of a third family of compact stars which may exist besides 
the two known families of white dwarfs and neutron stars. We show 
how an experimental proof of the existence of a third family by mass 
and radius measurements may provide a unique signature for a phase 
transition inside neutron stars. 

PACS: 12.38.Mh, 26.60,+c 

Keywords: Quark matter; Neutron stars; Third family 



1 Introduction 

The properties of nuclear matter under extreme conditions are subject to in- 
tense experimental and theoretical studies in heavy-ion physics as well as in 
nuclear astrophysics. Especially neutron stars provide a unique astrophysical 
environment to study the properties of matter in the region of low temper- 
atures and high baryonic densities. The densities deep inside neutron stars 
might be sufficient to populate numerous new forms of matter which may 
compete with each other. It is expected that at about two times normal nu- 
clear matter density hyperons appear in neutron star matter |], |2], ]3|, f|]. Also 
the formation of Bose-Einstein condensates |5], |6], |7j or the phase transition to 
deconfined quark matter have been considered |], |^, |], H 0, |ll], [12"| . For a 



recent review on possible phase transitions in neutron stars see |13| . However, 
our understanding of the properties of neutron star matter is still restricted 
due to theoretical limits and uncertainties in modeling the equation of state 
(EOS). Measurements of neutron star masses and radii may put important 
constraints on the EOS at large density. The masses of about twenty neutron 
stars (mainly pulsars in binary systems) have been measured and found to 
be quite near the "canonical" mass of M = 1.4 M & (M is the solar mass) 



14}| . Unfortunately, the currently known masses do not constrain the EOS 
unambiguously. While some soft EOS - which can only account for smaller 
masses - can be ruled out, it is not possible to differentiate between the class 
of stiffer EOS which give higher neutron star masses. However, these EOS 
can be classified by their different predictions for neutron star radii. Precise 
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measurements of neutron star radii can therefore give important information 
on the property of matter and possible exotic phases at large density. Some 
indirect radius estimates from observations of neutron stars seem to indicate 
surprisingly small radii. Recent radius estimates for Geminga (PSR 0630+17) 



by Golden and Shearer |]T5[ suggest an upper bound of the apparent radius of 
Roo < 9.5]t|; km for an assumed blackbody source and < lO.Ol^km with 
the presence of a magnetized H atmosphere. Walter and Matthews obtained 
i? QO <14km for the recently discovered isolated neutron star RXJ1856. 5-3745 
[ |i~6H . Generally the apparent radius -Roo at a distance d — > oo is related to the 
local stellar radius R by 



= R/yJl - 2GM/Rc 2 . (1) 

If we assume a mass of M = 1.4 M Q an apparent radius R^ of about 11 — 13 km 
would restrict the radius to the quite low values of R ~ 7 — 10 km. (Due to 
(HD, radii below i?^ ~ 11km would imply masses below 1.4 M .) A radius 
of i?<10km for a typical neutron star mass is furthermore consistent with 



the findings by Haberl and Titarchuk [17| for the X-ray burster 4U 1820-30 



with the findings by Pavlov and Zavlin [H| for the millisecond pulsar J0437- 



4715 and by Li et al. [pl| for Her X-l. (In the context of Ref. |19j see also 



Ref. |2C| and Refs. pj, If radius estimates of -R^lOkm indeed should 

prove to be reliable in connection with typical neutron star masses, the care- 
ful interpretation of these results in terms of neutron star matter properties 
would be an intriguing challenge for nuclear astrophysics. Since the appear- 
ance of a deconfined quark matter phase inside neutron stars is one possible 
mechanism to soften the EOS (and therefore to reduce the radius of the star), 
we want systematically address the question to which amount the compact- 
ness of a neutron star can be attributed to the presence of a quark phase. 
Such investigations are necessary also in view of more "exotic" interpretations 



of small radii e.g. in terms of "strange stars" |L9| , |23| , gj}] . It is important 
to study whether such interpretations are a consequence of ruling out other 
mechanisms suitable to soften the EOS at high densities (e.g. by hyperons, 
kaon condensates or quark phases) or provide only one possible explanation of 
the observed neutron star properties. For this purpose we will apply various 
hadronic EOS in the relativistic mean-field model (including hyperons) and a 
wide range of model parameters of the quark phase which we will describe in 
an extended MIT bag model. 

The EOS of the hadronic and the quark phase will be discussed in Sec. |2|. 
In Sec. [3] we discuss the construction of a first-order phase transition from 
hadronic to quark matter and present the results for the EOS. In Sec. f| we 
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apply the EOS to the calculation of the neutron star structure and discuss the 
influence of quark phases. In Sec. [| we discuss the properties of a possible 
"third family" of compact stars. Such hypothetical family of compact stars 
might exist besides the two known families of white dwarfs and neutron stars 



25, Gq, ETA. In this work we show how future measurements of masses and 



radii of only two compact stars may prove the existence of a third family. 
Moreover, we argue that this would provide a novel signature for a phase 
transition inside neutron stars, revealing a particular behavior of the EOS in 
the almost unknown regime of supernuclear densities. Finally, we summarize 
our findings in Sec. 0. 



2 Equations of state 

We want to start with a discussion of the different models used to describe 
the equation of state of "neutron star matter" i.e. cold, charge neutral mat- 
ter in beta equilibrium. Neutron star matter covers a wide range of densi- 
ties. Starting from the density of iron (e pa 8g/cm 3 ) which builds the surface 
of a neutron star, densities up to few times normal nuclear matter density 
(eo = 140MeV/fm 3 pa 2.5 x 10 g/cm 3 ) can be achieved in the center of the 
star. Unfortunately there is yet no single theory capable to meet the demands 
of the various degrees of freedom opened up in neutron star matter in different 
densities regimes. We are therefore forced to use different models at different 
density ranges. 



2.1 Subnuclear densities 



For subnuclear densities we use the Baym-Pethick-Sutherland EOS |28|, which 
describes the crust of the neutron star. Up to densities of e ^10 7 g/cm 3 matter 
is composed of a Coulomb lattice of 2gFe nuclei. The pressure is dominated 
by degenerate electrons. At higher densities nuclei become more and more 
neutron rich and at e^rip ~ 4 x 10 11 g/cm 3 the most weakly bound neutrons 
start to drip out of the nuclei. For a detailed discussion of the Baym-Pethick- 
Sutherland EOS see | 29| . A review of the properties of neutron rich matter at 

m. 



subnuclear densities can be found in 



2.2 Nuclear densities 

The Baym-Pethick-Sutherland EOS is matched at e ~ 10 14 g/cm 3 pa eo to a 
hadronic EOS calculated in the framework of the relativistic mean-field (RMF) 
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model. At these densities nuclei begin to dissolve and nucleons become the 
relevant degrees of freedom in this phase. The RMF model is widely used 



for the description of dense nuclear matter [31, 32, RL especially in neutron 



stars. For an introduction into the RMF model see e.g. 0. We use three 
EOS calculated by Schaffner and Mishustin in the extended RMF model || 
(i.e. TM1, TM2, GL85) and one by Ghosh, Phatak and Sahu @. For the 
latter one we use GPS as an abbreviation. These models include hyperonic 
degrees of freedom which typically appear at e « 2 — 3e . Table [I] shows the 
nuclear matter properties and the particle composition of the four EOS. Even 



Hadronic EOS 


TM1 


TM2 


GL85 


GPS 


reference 


i 


I 


@ 


H 


Po [fm" 3 ] 


0.145 


0.132 


0.145 


0.150 


B/A [MeV] 


-16.3 


-16.2 


-15.95 


-16.0 


K [MeV] 


281 


344 


285 


300 


m* N /m n 


0.634 


0.571 


0.770 


0.830 


a sym [MeV] 


36.9 


35.8 


36.8 


32.5 


composition 


a 


a 


a 


b 


a) n, p, e", fu~ 


, A, S-, £ 


v+ 







b) n, p, e ,n , A, £ 



Table 1: Nuclear matter properties of the hadronic EOS. The saturation den- 
sity and the binding energy is denoted by po an d B/A, the incompressibility 
by K, the effective mass by m* N /m n and the symmetry energy by a sym . The 
particle compositions are shown at the bottom of the table. 

if the relevant degrees of freedom are specified (in our case basically nucleons 
and hyperons) the high density range of the EOS is still not well understood. 
The use of different hadronic models should reflect this uncertainty to some 
degree. In the following we will denote the phase described by the Baym- 
Pethick-Sutherland EOS and by the RMF model as the hadronic phase (HP) 
of the neutron star. 

2.3 Quark matter 

At densities above eo we allow the HP to undergo a phase transition to a 
deconfined quark matter phase (QP). This section describes how we model the 
QP in weak equilibrium. The QP consists of u, d, s quarks and electrons in 
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weak equilibrium i.e. the weak reactions 

d — > u-\-e~ + v e -, (2) 
s — > u + e~ + v e - , 
s + u < — > d + u, 

imply relations between the four chemical potentials p u , pa, A*s, Pe which read 

Ps = Pd = Pu + Pe- (3) 

Since the neutrinos can diffuse out of the star their chemical potentials are 
taken to be zero. The number of chemical potentials necessary for the descrip- 
tion of the QP in weak equilibrium (the number of components) is therefore 
reduced to two independent ones. We choose the pair (p n , p e ) with the neutron 
chemical potential 

Pn = Pu + %Pd- 

In a pure QP (in contrast to a QP in a mixed phase, which we will discuss later) 
we have to require the QP to be charge neutral. This gives us an additional 
constraint on the chemical potentials 

PT = \pu - ^Pd ~ ^Ps ~ Pe = 0. (4) 

Here p® p denotes the charge density of the QP and pf (/ G u,d,s), p e the 
particle densities of the quarks and the electrons, respectively. The EOS can 
now be parametrized by only one chemical potential, say p n . At this point 
it should be noted that the arguments given here for the QP also holds for 
the HP. There one also ends up with two independent chemical potentials 
(e.g. p n , p e ) if one only requires weak equilibrium between the constituents of 
the HP, and with one chemical potential (e.g. p n ) if one additionally requires 
charge neutrality. As we will discuss later, the number of independent chemical 
potentials plays a crucial role in the formulation of the Gibbs condition for 
chemical and mechanical equilibrium between the HP and the QP. 

To calculate the EOS of the QP we apply the "effective mass bag model" 



11], Q . In this model, medium effects are taken into account in the framework 
of the MIT bag model |35], by introducing density-dependent effective quark 
masses. The effective quark masses follow from the poles of the resummed 
one-loop quark propagator at finite chemical potential (see Fig. |l|). Details of 
the calculation of the effective quark masses in the so called hard dense loop 



approximation can be found in [M. One finds [|34]. 37 



m 



m f imj g 2 p} 



4 6vr 
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with the current quark mass m/ of the quark flavor / and the strong coupling 
constant g. These effective quark masses show the behavior of increasing mass 
with increasing density (chemical potential). For the current quark masses we 
choose m„ = irtd = and m s = 150 MeV. The effective masses @ can now be 
used in the Fermi gas expression for the particle density 



3/2 



PfM = ^2 K~ m / M\ > ( 6 ) 

where df is the degree of degeneracy |38|. The pressure pj and energy density 
€f can be extracted from the thermodynamic relations 

= dPf d ^ and e f (n f )+pf(iif) = (?) 

If g —>■ the effective quark masses © approach the current quark mass rrif. 
In this limit, the EOS equals to a free MIT bag Fermi gas EOS. We therefore 
can use g to study the deviations from the free Fermi gas owing to the influence 
of medium effects. For this purpose we take g as a parameter ranging from 



g = (no medium effects) to g ~ 2 — 3 (a s ~ 0.31 — 0.72). (It was found in [11 
that a phase transition to a pure QP does not occur inside a typical neutron 
star for g^2. Larger values are therefore of minor interest for our purpose 
of studying the influence of quark phases on the properties of neutron stars.) 
The overall energy density cqp and pressure pqp of the QP is the sum over 
all flavors plus the Fermi gas contribution e e , p e of the uniform background of 
electrons plus the contribution of the phenomenological bag constant B. 

£qp = e u + e d + e s + e e + B, (8) 
Pqp = Pu+Pd + Ps+Pe- B, 

with 

The bag constant B is introduced in the usual way by adding it to the energy 
density and subtracting it from the pressure. Due to the phenomenological 
nature of the bag model, an exact value of B is not known. Therefore we will 



• > = — > — + 

Figure 1: Resummed quark propagator. 
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take B free parameter ranging from B 1 / 4 = 165 MeV (B « 96 MeV/fm 3 ) 
to 5 1/4 « 200 MeV (B « 208 MeV/fm 3 ). The lower bound of the bag constant 
is reasonable if one requires that the deconfmement phase transition to the 
QP should not occur at densities below e m €q. We do not consider here 
the possibility of the existence of pure QP stars, so called "strange stars" 
or the related case of an absolutely stable QP (stable strange quark matter) 
with respect to the HP as it was supposed by Witten |39|. These possibilities 
are in detail discussed in the literature and we refer to |T|, |2], ^0], f|T| for an 
overview. To find the QP to be absolutely stable, their energy per baryon E/A 
in equilibrium (pqp = 0) must be lower than the energy per baryon of 56 Fe 
which is about 930 MeV. This requires small bag constants of B 1 / 4 ^ 155 MeV 
(if we neglect medium effects). At our lower bound of the bag constant of 
B 1 / 4 = 165 MeV we obtain an energy per baryon of E/A m 990 MeV which 
is about 60 MeV unbound with respect to 56 Fe. As one can imagine from (§), 
the bag constant plays a crucial role in the question whether or not a QP can 
exist in the interior of neutron stars. 



3 Phase transition 



Using the four EOS of the HP discussed in section |2J] and |272| and the EOS of 
the QP discussed in the last section ( |2.3|) , we study the deconfmement phase 
transition from the HP to the QP and calculate the corresponding total EOS. 
In particular we are interested in the influence of the model parameters of the 
QP (i.e. the bag constant B and coupling constant g) on the phase transition 
densities. (The existence of a QP inside the neutron star clearly requires the 
transition density to be smaller than the central density of the star.) Before 
we come to this, we have first to discuss how to construct a first-order phase 
transition from the confined HP to the deconfined QP. 

3.1 Construction of the phase transition 

We calculate the phase transition according to Glendenning who has first re- 
alized the possibility of the occurrence of a mixed phase (MP) of hadronic and 
quark matter in a finite density range inside neutron stars |1], Q . For simplic- 



ity we neglect Coulomb and surface effects in the MP as studied in Ref. [}43 
The essential point in the calculation of the MP is that total charge neutrality 
can be achieved by a positively charged amount of hadronic matter and a neg- 
atively charged amount of quark matter. As already discussed in Sec. ^]3| we 
have to deal with two independent chemical potentials (// n , /x e ) if we impose 
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the condition of weak equilibrium. Such a system is called a two- component 
system. The Gibbs condition for mechanical and chemical equilibrium at zero 
temperature between the both phases reads 

PHp(l^n,l^e) = PQp{Hn, He)- (10) 

Using Eq. fliCf ) we can calculate the equilibrium chemical potentials of the 
MP where pup = Pqp = Pmp holds. Fig. ^| illustrates this calculation for a 
specific choice of the model parameters (see figure caption). The HP— >MP 




Figure 2: Gibbs phase construction of a two-component system. Plotted is the 
pressure of the hadronic phase pnp and the pressure of the quark phase pqp 
as a function of the two independent chemical potentials fi e . The white 
lines HP and QP on the pressure surfaces show the pressure of the hadronic 
phase and the quark phase under the condition of charge neutrality. At low 
pressure matter is in its charge neutral HP. At point CD the HP reaches the 
intersection curve MP of the mixed phase. This curve is the solution of the 
Gibbs condition (|T0"D. At a pressure above point [2] matter consists of a pure 
QP. EOS of the HP is GPS, EOS of the QP uses S 1 / 4 = 170 MeV and g = 2. 

phase transition takes place if the pressure of the charge neutral HP (white 
line) meets the pressure surface of the QP. This happens at point Q] in Fig. ||. 
Up to this point the pressure of the QP is below the pressure of the HP, 
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making the HP the physically realized one. Up to point [2] the phase follows 
the MP curve which is given by the Gibbs condition (0). At point [2] the MP 
curve meets the charge neutral QP curve (white line) and the pressure of the 
QP is above the pressure of the HP, making the QP the physically realized 
one. The resulting chemical potentials are again shown in Fig. |3]. Below the 
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Figure 3: Electron chemical potential function of the neutron chemical 

potential p n . Up to /j$ matter is in its charge neutral HP. From //W to /ij^ 
matter consists of a mixed phase (MP) of a charge positive HP (p^ p > 0, p c 
denotes the charge density) and a charge negative QP (pf p < 0). The MP 
curve is calculated from the Gibbs condition (110) . Above matter consists 
of a pure charge neutral QP. EOS of the HP is GPS, EOS of the QP uses 
B 1 ^ = 170 MeV and g = 2. 

charge neutral HP curve and above the charge neutral QP curve the HP is 
positively charged (p^ p > 0, p c denotes the charge density) and the QP is 
negatively charged {p® p < 0). Therefore, the charge of hadronic matter can 
be neutralized in the MP by an appropriate amount of quark matter. Up 
to the critical chemical potential p$ the electron chemical potential of the 
charge neutral HP is increasing. (This is mainly because electrons have to 
neutralize the increasing amount of protons appearing in the HP.) At /iW the 
phase transition to the MP takes place and p e can be reduced due to the 
first appearance of negatively charged quark phase droplets immersed in the 
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hadronic phase. (For a detailed discussion of the geometrical structure of the 
MP see U . ) For every point on the MP curve one has to calculate the volume 
proportion 

VQP + VHP 

occupied by the quark phase by imposing the condition of global charge neu- 
trality in the MP 

XP? P + (1-X)P? P = 0. (12) 
From this, the energy density €mp of the MP follows as 

^-MP = X^QP + (1 - x) e HP- (13) 

In the MP the volume proportion of the quark phase is monotonically increas- 
ing from x — at /j$ to x — 1 & t /4i where the transition into the pure QP 
takes place. Taking the charge neutral EOS of the HP (Sec. |2.2|) for p, < //W, 
Eq. (ffg), (TT2]) and ([□]) for the MP for //W < p n < /4 2] and the charge neutral 
EOS of the QP (Sec. |2.3|) for //„ > p$ we can construct the full EOS in the 
form p = p(e). For simplicity we denote the complete EOS as the hybrid star 
EOS. 

As one realizes from Fig. ^ the phase transition construction in a two- 
component system leads to a monotonically increasing pressure of the MP 
with increasing density [p2]| . This is quite different from the constant pressure 
in a one-component system (with only one independent chemical potential) 
which neutron star matter is occasionally assumed to be. For the existence 
of a MP inside the star an increasing pressure is crucial. In turn, such a MP 
allows for a softer EOS over a wide density range which would be excluded in 
the one-component treatment. We will discuss this point in more detail in the 
appendix. There we show that a simplified treatment of the phase transition 
with only one independent component results in an overestimation of the phase 
transition density. 



3.2 Results for the equation of state 

In this section we present the results of the phase transition construction dis- 
cussed in the last section. We study in particular the influence of different 
hadronic EOS and the influence of the model parameters of the QP on the 
properties of the phase transition. Fig. |] shows the low density range of the 
hybrid star EOS in the form p = p(e) for the four hadronic EOS and two cou- 
pling constants of the QP. The bag constant is chosen to be B 1 / 4 = 170 MeV. 
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Figure 4: Low density range of the EOS for several hadronic EOS and two 
coupling constants g of the quark phase. The bag constant of the quark phase 
is chosen as B lj/4 = 170 MeV. Shaded regions correspond to the mixed phase 
(MP) parts of the EOS. HP denotes the hadronic phase, eo = 140 MeV/fm 3 . 



The MP part of the EOS is shaded gray. As already noted, even in the MP the 
pressure is monotonically increasing with increasing density. For that reason 
the MP does exist in the interior of the neutron star if its central energy den- 
sity exceeds the HP^MP transition density. One can see that the HP— >MP 
transition density depends only slightly on the choice of the hadronic EOS. 
This reflects the rather small uncertainty in the hadronic EOS in the den- 
sity range up to about 1.5 eo where hyperonic degrees of freedom are still not 
present. (In the HP the hyperons typically appear at 2—3 e [0.) The HP— >MP 
transition density also depends only weakly on the influence of medium effects 
which shift the transition densities to slightly higher values. This was already 
found in |IT[] . Obviously the EOS gets softer compared to the HP due to the 



onset of the phase transition. (By "softer" we denote a smaller pressure at 
a fixed energy density e). In Fig. [| we show the same hybrid star EOS in 
the high density region. In contrast to the HP— >MP transition density we 
see that the MP— >QP transition density is quite sensitive to the choice of the 
hadronic EOS as well as to the influence of medium effects in the QP. As- 
suming e.g. a star with a central energy density of e = 5 eo we see that no 
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Figure 5: High density range of the hybrid star EOS of Fig. |j. The bag 
constant of the quark phase (QP) is chosen as B 1 ^ = 170 MeV. Shaded regions 
correspond to the mixed phase (MP) parts of the EOS. eo = 140MeV/fm 3 . 



conclusive statement can be made about the composition of the center of the 
star. Neglecting medium effects (g = 0), TM1 and TM2 predict a QP core 
(QC) while GPS and GL85 predict a MP core (MC). Taking medium effects 
into account (g = 2) only TM2 predicts a QC. The strong sensitivity of the 
transition densities can qualitatively be understood if we look again at Fig. [2|. 
In the high pressure region where the MP— >QP transition takes place, the HP 
and the QP pressure surfaces are nearly parallel. Therefore a small change in 
one EOS is able to produce considerably large changes in the MP— >QP tran- 
sition density. Furthermore one should note in Fig. |5| that (despite a change 
in the phase transition densities) the influence of medium effects on the pure 



QP EOS in the form p = p{e) is negligible. This was already found in |34 



Up to now we have neglected the influence of the bag constant B which 
enters in the QP EOS. Since B enters negatively into the QP pressure in Eq. 
(H) a larger B will shift the pressure surface PQp(fJ, n , /-0 m Fig- @ down which 
moves the MP intersection curve to higher pressures. This finally also shifts 
the transition densities to higher values. Fig. |6| and Fig. [7| demonstrate this 
behavior for different hybrid star EOS using GPS and TM2 for the HP (the 
other hybrid star EOS using TM1 and GL85 show a very similar behavior). 
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Figure 6: Phase transition densities for three different coupling constants g 
(thin lines) as a function of the bag constant B. The shaded region (MP) 
marks the density range of the MP at g = 0. The thick lines show the central 
energy density of a 1.4M star. As long as the thick line is to the right of the 
corresponding MP range (e.g. the shaded region for g = 0) the star possesses 
a pure quark phase core, eo = 140MeV/fm 3 . 

In these figures the gray shaded area corresponds to the density region of 
the MP for g = (i.e. no medium effects). Obviously the HP— >MP and the 
MP— >QP transition densities are increasing with increasing B and increasing 
g where the MP— >QP transition density again is more sensitive. An increase 
of B l l A by only 5 MeV is able to increase the MP-^QP transition density 
by about one eo- From Fig. ^| it is furthermore interesting to note that the 
HP^MP transition density seems to be only slightly sensitive to a change of 
the model parameters as long as its density is below e ~ 2e and strongly 
sensitive above (comparable to the MP— >QP transition density). This effect 
seems to be related to the appearance of hyperons (mainly S~ and A) in the 
HP EOS at densities of e « 2 — 3eo which has already been noted in ||11| . 

It is important to note that the increase of the transition densities with 
increasing B does not necessarily disfavor the existence of deconfined phases 
inside the star since also the central energy density can be increased by chang- 
ing this parameter. This is illustrated by the thick lines in Fig. ^| and Fig. [7] 
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Figure 7: Phase transition densities for three different coupling constants g 
(thin lines) as a function of the bag constant B. The shaded region (MP) 
marks the density range of the MP at g — 0. The thick lines show the central 
energy density of a 1.3M star, eo = 140MeV/fm 3 . 

which show the central energy density of a star at fixed mass as a function 
of the bag constant. (The calculation of neutron star masses is discussed in 
the next section. In Fig. we choose a mass of M = 1.3 M instead of the 
"canonical" value of M = 1.4 M & since it turns out that for some (B, g) pa- 
rameters the maximum mass (M max ) is slightly below 1.4 M Q .) At least for 
g = and g = 1 the central energy density is increasing with B for low values. 
From these figures we can already see which kind of phase (HP, MP or QP) 
is realized in the center of the star. E.g. from Fig. || we find that at g = 
the central energy density (thick solid line) of a M = 1.4M star is above the 
MP— >QP transition density for bag constants below B 1 / 4 « 180 MeV. In this 
parameter range the neutron star possesses a pure QP core (QC). At higher 
B up to B 1 / 4 « 195 MeV the star possesses a MP core (MC). At still higher 
B the star is purely hadronic (HC). We will come back to a discussion of the 
(B, g) parameter range leading to one of the possibilities (QC, MC, HC) in 
the next section. 

In the following we first want to focus in some more detail on the interesting 
behavior of the central energy density with increasing B. The central energy 



15 



density of a star at fixed mass depends on the softness or stiffness of the 
underlying EOS. Soft EOS possess a larger central energy density as compared 
to stiff EOS. In a simple way one could imagine that stars corresponding to 
soft EOS are compressed to smaller radii and therefore possess a larger central 
density. In our case (e.g. Fig. || thick line g = 0) two competing effects are 
important. By increasing B, i) the low density range of the hybrid star EOS 
gets stiffer due to an enlarged density range of the HP (which is stiffer than 
the MP) while ii) the density range of the QP gets softened directly by a larger 
bag constant. Effect ii) dominates for smaller B since large parts of the star 
are made of the QP. Therefore the central energy density is increasing with 
B. In both figures this holds up to about B 1 ^ m 175 MeV for g = 0. We will 
see in the next section that in this range the internal structure of the star is 
only slightly sensitive to a change of B because both the MP— >QP transition 
density and the central energy density are equally increasing with B. The 
situation is reversed for bigger B and the central energy density is decreasing 
towards the value of the pure hadronic star. This behavior demonstrates two 
things. Firstly, the central energy density depends quite sensitively on the 
softness or stiffness of the EOS which effect the composition of neutron star 
matter. Thus it is not possible to talk about a "typical central energy density" 
of a neutron star without having a particular EOS in mind. Secondly, the 
increase of the bag constant - which makes the QP energetically less favorable 
- does not automatically disfavor the existence of a QP inside the star since 
also the central energy density is increased. Indeed, there is a region of the 
bag constant (_B 1//4 £175 MeV) where the internal structure of the star depends 
only slightly on it. We will come back to this point below. We furthermore 
should note that the discontinuous fall of the central energy density in Fig. [7] 
for g = at B 1 ^ « 180 MeV and for g — 1 at slightly lower B is related to 
the existence of "neutron star twins" which we will discuss in Sec. |5]. 



4 Neutron star structure 

With the evaluated hybrid star EOS presented in the last section we now turn 
to analyse the structure of the corresponding non-rotating neutron stars by 



solving the Tolman-Oppenheimer-Volkoff (TOV) equations [44 . These equa- 



tions describe the balance between the gravitational force and the pressure 
given by the Fermi pressure of the particular EOS. This leads to a relation be- 
tween the mass M and the radius R of the neutron star in general relativistic 
hydrostatic equilibrium. 
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Figure 8: Mass-radius relations using the GPS EOS for the HP. The uppermost 
line denoted by H corresponds to the pure hadronic solution. All other curves 
allow for a phase transition to the quark phase. The small numbers denote 
the coupling constant g employed for the QP (dashed lines denote steps of 0.5 
in g). Objects in the gray shaded MC region posses a mixed phase core while 
stars in the darker shaded QC region posses a pure quark phase core. Note 
the radius separation of about 2 — 4 km between the stars possessing a quark 
core and the pure hadronic stars of the same mass. 



4.1 Mass-radius relations 

Fig. § and Fig. |9] show the mass-radius (MR) relations for two hybrid star 
EOS (GPS and GL85) and several bag constants and coupling constants of 
the QP. The MR-relation denoted by H corresponds to the solution for the 
pure hadronic star (without any phase transition). The central energy density 
e c of the stars is increasing from large radii and small masses (down to the 
right) to the upper left edge where the critical central energy density e c ™ 1 is 
reached at the maximum mass and the minimal radius of the star. From the 
two differently shaded regions denoted by QC (QP core) and MC (MP core) 
one sees whether e c is large enough to possess a QC or a MC in the center of the 
star. For small bag constants (-B 1//4 <180 MeV) and small coupling constants 
(lower curves) e c is sufficient to exceed the MP— >QP transition density before 
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Figure 9: Mass-radius relations using the GL85 EOS for the HP. The upper- 
most line denoted by H corresponds to the pure hadronic solution. All other 
curves allow for a phase transition to the quark phase. The small numbers de- 
note the coupling constant g employed for the QP (dashed lines denote steps 
of 0.5 in g). Stars in the gray shaded MC region posses a mixed phase core 
while stars in the darker shaded QC region posses a pure quark phase core. 



e^. rit is reached (cf. Fig. ||). This leads to a pure QC inside the star. For a 
more detailed discussion of MR-relations and in particular of the influence of 



medium effects on the MR-relation see [|TT| . 

It is now our major aim to see what the effect of a pure QC on the mass 
and the radius of a neutron star is. Concerning the mass of the star we find 
that M max of a pure hadronic star (H) is reduced by about 0.1—0.3 M due to 
the phase transition and the corresponding softening of the EOS. Nevertheless, 
both hybrid star EOS are able to explain a typical neutron star mass of about 
1.4 Mq. The most compact objects seem to have radii of about 9 — 10 km. In 
our model calculations these radii are only reached in stars possessing a pure 
QC (especially for g m 0). Very similar results are found for the two other 
hybrid star EOS (TM1 and TM2) which are not shown here. If we furthermore 
compare the radii of stars located in the QC region with the radii of the pure 
hadronic stars (H) of the same mass we can see that the QC stars are about 
2— 4 km smaller than the corresponding hadronic stars. Hence, a neutron star 
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Figure 10: Compactness ([14]) of neutron stars possessing a pure quark core with 
respect to the pure hadronic star of the same mass. Obviously the presence 
of quark phases inside the star reduce the radius by typically 20 — 30%. The 
shaded areas stem from varying the model parameters of the quark phase (bag 
constant B (shading) and coupling constant g) which lead to a pure quark 
core inside the star. Especially for GPS, GL85 and TM1 the compactness 
ratio seems to be rather insensitive to the particular hadronic EOS and to the 
model parameters of the quark phase. 



is about 20—30% more compact due to the presence of a pure QP in its center. 
We define the compactness of a neutron star possessing a QC as 

( ^ _ R H {M) - Rq C {M) 
compactness(M) = — ^ , (14) 

with the radius Rqc of the QC star and the radius Rh of the pure hadronic 
star of the same mass (cf. Figs. || |9|). The compactness therefore corresponds 
to the ratio to which the radius of a pure hadronic star can be reduced due to 
the presence of a pure QP. The shaded areas in Fig. |1^ come from calculating 
the compactness for all stars located in the QC region of the corresponding 
MR-relation (cf. Figs. §, ||) for different bag constants (which are depicted 
in Fig. [iy by different shadings). It is remarkable that the compactness of 
0.2 — 0.3 depend only slightly on the particular HP EOS and on the model 
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parameters applied in the QP. Especially an increase of the bag constant B 
(darker shading) only narrows the mass range leading to a QC. At the same 
time the range of the compactness is only slightly altered. A reduction of the 
radius of an hadronic star by 20—30% due to the presence of a pure QP seems to 
be a typical value. Only for g « (which in general leads to the most compact 
stars) and the TM2 hybrid star EOS a radius reduction of about 40% seems 
possible (cf. Fig. [H|). This might be related to a rather large incompressibility 
K, a small saturation density po> an d a small effective mass m* N of the TM2 
EOS compared to the other EOS (cf. Tab. []]). 

4.2 Internal structure 

In this section we further study how the internal structure of the star e.g. the 
radius of the quark core (QC) or the thickness of the MP depends on the model 
parameters of the QP. Fig. [TT] and Fig. [12] show the neutron star cross sections 
for a range of bag constants B and a fixed neutron star mass. The respective 
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Figure 11: Internal structure of a M = 1.4M Q neutron star as a function of 
the bag constant. Left panel g — (no medium effects), right panel g = 1.5. 
The hadronic EOS is GPS. 

left panels show the results without medium effects (g = 0) while the right 
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Figure 12: Internal structure of a M = 1.3M neutron star as a function of 
the bag constant. Left panel g = (no medium effects), right panel g = 2. 
The hadronic EOS is TM2. 

panels include medium effects in the QP. The three differently shaded regions 
QP, MP and HP again correspond to the parts of the star made of quark phase, 
mixed phase and hadronic phase, respectively. The thick line to the right of 
the HP marks the surface of the star. It is obvious that basically only the 
presence of a pure QP is able to reduce the radius of the star significantly. For 
g = (left panel) we can see in both figures that only for 5 1//4 <180 MeV a 
pure QC can exist. For B l ^< 175 MeV the radius of the QC can be as large 
as 4 — 7 km at an overall radius of about 10 km. For these values of B we can 
see that the internal structure of the star (especially the overall radius) only 
slightly depends on a change of B. This behavior can be traced back to the 
findings of Sec. |3.2| . There we found that despite of an increase of B from 
B 1 ^ = 165 MeV to 175 MeV (which increases the MP^QP transition density 
in the GPS case by about 3 e ) a pure QP is not excluded due to a similar 
increase of the central energy density e c (cf. Figs. |] and [7]). 

On the right panels we can see how the influence of medium effects change 
the internal structure of the star. Since the HP— >MP and MP— >QP transition 
densities at fixed B increase with increasing g (cf. Fig. ^|and [?D medium effects 
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Figure 13: Parameter range (bag constant and coupling constant) of the in- 
ternal structure of a M = 1.4M neutron star. Parameters in the QC region 
lead to neutron stars with pure quark core while MC or HC corresponds to a 
mixed core or a pure hadronic star, respectively. The dashed lines are lines 
of constant radius given as the percentage of the radius of the hadronic star. 
The hadronic EOS is GPS. 



are able to transform a star with MC into an pure hadronic star (e.g. Fig. [11 



B 1 ' 4 = 190 MeV) or a star with QC to one with a MC (e.g. Fig. 0, B 1 ^ 
175 MeV). Only a lower bag constant can compensate for this effect to a certain 
degree. The solutions for the TM1 and GL85 hybrid star EOS (which are not 
shown here) lead to very similar results. 

Finally, we illustrate the complete (B, g) parameter range in Fig. [13] and 
|l4l to see which values lead to which kind of internal structure (QC, MC or 
HC (pure hadronic star)). Again we can see that - if we want to keep the 
internal structure of the star nearly fixed - the influence of medium effects can 
only be compensated by lowering the bag constant. The dashed lines show the 
neutron star radii in terms of the radius of the pure hadronic star. These lines 
correspond to lines of constant compactness (0). In both figures the borderline 
between the QC and the MC region appears at about 80% (corresponding to a 
compactness of 0.2). This demonstrates once more that - almost independently 
of the model parameters of the QP - a neutron star with a MC can at most be 
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Figure 14: Same as Fig. \LQ but for a neutron star mass of M = 1.3M . The 
hadronic EOS is TM2. 



about 20% more compact than a hadronic star of the same mass while a star 
with a QC is typically 20 — 30% more compact than a hadronic star. 
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5 Phase transitions and the third family 



5.1 Gerlachs criterion 



Almost 30 years ago, it was found by Gerlach |26| that a third family of 
stable equilibrium configurations of compact stars are not forbidden by general 
relativity. Such a third family could therefore in principle exist besides the two 
known families of white dwarfs and neutron stars. (See Fig. ^ for a schematic 
view of the three families in a MR-radius relation. The criteria for stability 
and instability of the shown ranges will be discussed in Sec. |5.3| .). Gerlach 
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Figure 15: Schematic mass-radius relation showing three stable families of 
compact stars. The letters A, B, . . . , I refer to critical points (turning points) 
where a vibrational mode changes stability [2T3, The stability (solid lines) 
or instability (dotted lines) of the three lowest-lying modes (n = 0,1,2) is 
depicted by the numbers. Higher modes are stable. See text for more details. 



found that a necessary condition for the theoretical possibility of this family 
is a sufficiently large discontinuity in the speed of sound (c 2 s = dp/de) of the 
corresponding EOSQ. As the density is increased, the speed of sound has to 
increase abruptly at a density exceeding the central energy density of the 

lr inis was proved by means of an inversion of the Tolman-Oppenheimer-Volkoff equations. 
In this context see also the work by Lindblom p6[. 
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neutron star of limiting mass (point C in Fig. |15|). Rather smooth EOS (as 
comprehensively studied by Wheeler et al. [f45] ) does not possess the possibility 
of a third family. Thus a verification of the existence of a third family in nature 
is equivalent to the verification of a phase transition in the EOS of matter at 
large densities. 



5.2 Quark phases and the third family 

From Gerlachs criterion it remains however unclear i) which physical mecha- 
nism could produce a discontinuity (a phase transition) sufficient to support 
the stability of a third family and ii) whether a formation process exists in 
nature to physically realize these objects. Concerning i) it was recently found 
by Glendenning and Kettner |27| that certain EOS that describe a first-order 
phase transition to deconfined quark matter (like the EOS also considered in 
this work) could indeed lead to a third family. Due to partially overlapping 
mass regions of the neutron star branch and the branch of the third family 
it is possible that non-identical stars of the same mass can exist. Such pairs 
are refered to as "neutron star twins" [p7| . In such kind of EOS the necessary 
discontinuity in the speed of sound is produced by the MP^QP transition of 
the deconfinement phase transition. As can be seen in our model from Fig. ||], 
especially the TM2 hybrid star EOS possesses a large discontinuity due to a 
rather low speed of sound at the high density end of the MP. In a particular 
parameter range of the bag constant we will see that the TM2 hybrid star EOS 
indeed allows for the existence of a third family^. Fig. [16] shows MR-relations 



using the TM2 hybrid star EOS. They are shown for six different bag constants 
£i/4 f rom to 183 MeV. We furthermore assume g = (no medium effects 
in the QP). The solutions for B 1 ' 4 = 175 MeV and B 1 ' 4 = 183 MeV do not lead 
to a third family and thus are similar to QC and MC solutions found in Figs. |8| 
and Increasing the central density beyond the critical point does not restore 
stability (corresponding to range C-H-I in Fig. [15]). This is not the case for bag 
constants in the small intermediate range 176 MeV < B 1 ' 4 < 182 MeV. There 
the mass-radius relations are splitted by an instable (dotted) range which cor- 



responds to C-D in Fig. [15] This enables the possibility of neutron star twins 
as shown e.g. for the low-density and high-density twin of mass M = 1.36M Q . 
While the low-density twin (on the neutron star branch) has a MC the high- 

2 For the other hybrid star EOS we have checked more than 200 MR-relations following 
from different combinations of B and g without finding any further solutions including a 
third family. Since the parameter range for a third family seems to be quite small we might 
however ignored some solutions. 
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Figure 16: Mass-radius relation using the TM2 EOS for the HP for different 
bag constants B and g = 0. For this EOS a third family only exists in the 
narrow parameter range of 176 MeV < 5 1 / 4 < 182 MeV. Obviously the stars of 
the third family possess approximately the same minimal radii as an "ordinary" 
hybrid star with pure QC as given by _B 1//4 = 175 MeV. 



density twin of the third family is more compact and possesses a QC. This is 
a general characteristic of neutron star twins PT| . The neutron star branch 
terminates in the MP owing to a small adiabatic index T = (e + p)/p ■ dp/de 
(cf. Fig. |^) P7| . Thus, all stars on the neutron star branch are MC stars (or 
HC stars at lower e c ). Only the larger adiabatic index of the QP can restore 
stability and therefore enables the QC stars of the third family. The internal 
structure of the neutron star twins with M = 1.36M is schematically shown 
in Fig. [17, 

In this section we have seen, that from the deconfmement phase transition 
described in our model a sufficiently large discontinuity in the speed of sound 
arises to allow a third family of compact stars. This possibility was first real- 



ized by Glendenning and Kettner ||27|| . Despite the open questions concerning 
the existence of a formation process, we now want to ask for the theoretical 
possibility of identifying a third family by means of mass and radius measure- 
ments of neutron stars. To do so we have first to discuss the stability criteria 
which ultimately separate the third family branch of a MR-relation (range D- 
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Figure 17: Schematic gross structure of the non identical neutron star twins 
marked in Fig. [TBI. Both stars have M = 1.36M and belong the same EOS 
(TM2 hybrid star EOS, B 1 ' 4 = 182 MeV, g = 0). 



E in Fig. [15]) from the neutron star branch (range B-C) in the same way as 
neutron stars are separated from white dwarfs. 



5.3 Stability criteria 

While hydrostatic stability of stellar configurations is assured by means of 



the Tolman-Oppenheimer-Volkoff equations [I3J a prove of dynamical stability 
requires an additional analysis of the radial vibration modes (acoustical modes) 
of the star. Due to dynamical instabilities no stable stars can have central 
densities in the range of about 10 9 — 10 14 g/cm 3 . This separates the family of 
the white dwarfs from the neutron star family. A configuration is stable against 
small perturbations around hydrostatic equilibrium if and only if the squared 
frequency eigenvalue uj^ of the fundamental vibrational mode is positive |2^, 
|45| . Dynamical stable and instable parts of a schematic MR-relation are shown 
in Fig. |15[ The letters A, B, . . . , I refer to critical points (turning points) where 
dM I de c = holds. The central energy density e c is rising from right to left 
along the curve A-C-I and A-C-G respectively. Reaching a critical point one 
radial vibrational mode (characterized by its number of nodes n) has to change 
stability. This means that the corresponding squared frequency eigenvalue 
changes its sign |29|, [|^] . (A stable mode requires a positive squared frequency 
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lo\ > 0). At critical points where dR/de c < holds (all points except H and F) 
an even mode changes stability while for dR/de c > (points H and F) an odd 
mode changes stability |29], fl5|] . Starting from a stable neutron star located in 
the range B-C (all lo\ > 0) we reach the critical point C by increasing e c . At 
point C an even mode changes stability due to dR/de c < and therefore gets 
instable. Since 

ul < u\ < u\ < - ■ ■ < u 2 n (15) 
holds ||, H this only can be the fundamental n = mode. A perturbation 



would cause such instable star to explode or collapse to a black hole. The only 
way to recover stability and therefore to make a third family possible at larger 
density is a further change of stability of the even n = mode at point D. (Note 
that point H can only change an odd mode). Due to ( |15"D no even mode expect 
n = can change stability at this critical point since all higher modes are bound 
by uj\ which is still positive. For that reason the region D-E corresponding to 
the third family is again stable against radial vibrations. Without a sufficiently 
large discontinuity in the speed of sound (i.e. for any reasonably smooth EOS) 
the mass-radius relation follows the curve C-H-I f45| , |47j . Increasing the central 
density then leads to a successive excitation of more and more unstable modes. 
(Only modes up to n = 2 are shown in Fig. [15|). 

As a important consequence of the stability analysis above, the mass and 
radius of point C has to be larger than the mass and radius of point D. There- 
fore it is possible that stars of the neutron star branch can in principle exist 
which are larger and heavier than stars of the third family. 

5.4 The third family as a possible signature for a phase 
transition 



Within our model we see from Fig. [16] that the masses and radii of the third 



family are comparable to the ones of the neig hbouring B 1 ^ = 175 MeV MR- 
relation which does not possess a third family. Moreover, for a radius and mass 
of R ~ 10 km and M « 1.35M Q we can see the interesting behavior that all 
branches of the third family intersect with the branch (B 1 ^ 4 = 175 MeV) of the 
neutron star family in just one point. Only the existence of such an example 
(even though found in a particular model) illustrates that the measurement 
of the mass and the radius of only one star (even if exact) could in general 
not provide the information to decide whether this star is a neutron star or an 
object of the third family. But how much points on a MR-relation do we in 
principle have to know (to measure) to decide whether or not a third family 
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R 

Figure 18: Schematic mass-radius relation showing a "rising twin" which can 
only exist if a third family of compact stars exists (AM <C AR). 

exists? The answer is two^\. Comparing two arbritay points on a neutron 
star MR-relation (e.g. £ 1/4 = 175 MeV) we see that the heavier star always 
possesses the smaller radius. But this must not hold for MR-relations including 
a third family as discussed in the previous section. Then it is possible that 
the heavier star possesses the larger radius. Such pair of stars always exists if 
a third family exists. For that it is necessary that the smaller star is located 
on the third family branch while the larger star is on the neutron star branch. 



This is schematically shown in Fig. 18 where we have used the same notation 



as in Fig. |15| We will refer to such pair as a rising twin. The idea of a rising 
twin as a signature for a third family is therefore based on the decrease of the 
mass with increasing radius of all MR-relations without a third family. This 
holds for almost all gravitational bound stars in which an increase of mass leads 
to a decrease of the radius due to the increasing gravitational attraction (see 
also the Figs. |9|). An exception of this is given by self-bound stars like the 
hypothetical strange stars which show the behavior of increasing radius with 
increasing mass (for small masses M oc R 3 holds) M. Also in some hadronic 
EOS (see e.g. p8| ) mass ranges of slightly increasing radius with increasing 
mass can exist. However, we can exclude these exceptions by additionally 
require that AM <C AR (G = c = 1) has to hold for a rising twin. While 
in the case of a third family even AM = is allowed (the neutron star twin) 
this is not possible for the exceptions for which typically AM ^> AR holds in 



3 For this it is important to note that the two masses and radii refer to two stars located 
on the same MR-relation. Therefore we can e.g. not refer to two stars with largely different 
temperature, magnetic field or rotational period. 



29 



the mass range of typical neutron stars. In the framework of our model (see 
Fig. [16]) the corresponding radius differences can grow as large as AR « 3 km 
at mass differences smaller than AM m 0.05 M e « 0.07km (M Q « 1.48km). 
For all possible rising twins of Fig. [16] AM is more than one order of magnitude 
smaller than AR. The requirement of small mass differences is furthermore not 
in contradiction to the experimentally known masses of neutron stars which 
are (at least for the classes of observed stars) in a remarkably narrow range of 
M = 1.35 ± O.O4M In particular the mass differences of the two stars in 



some double neutron star binaries seems to be quite smallQ 

In applying the stability considerations of Sec. |0] we conclude that the 
detection of a rising twin (by the measurement of masses and radii of two 
neutron stars) can prove the existence of a third family of compact stars. 
Gerlachs criterion (Sec. |5.1| ) shows that this is equivalent to the existence 
of a phase transition in the EOS at finite density. The deconfinement phase 
transition (as discussed in this work) is one possible scenario that might explain 
the hypothetical third family of compact stars (Sec. |5.2|). 



6 Summary and conclusion 

We have studied the properties of non-rotating cold neutron stars including 
the possibility of a phase transition to a deconfined quark phase (QP). To 
describe the confined hadronic phase (HP) of the star we have applied several 
hadronic EOS in the framework of the relativistic mean-field model (cf. Sec. |2.2| , 
|3"3]]). The QP was modeled by an extended MIT bag model (effective mass 
bag model) which includes medium effects by means of effective medium de- 
pendent quark masses (cf. Sec. pT3| , |34|, |i~l|). The influence of medium effects 
was parametrized by the strong coupling constant g. The construction of a 
first-order phase transition from confined to deconfined matter was performed 
by taking into account two independent chemical potentials (components) of 
neutron star matter (cf. Sec. [|, Fig. 0). This construction allows for the exis- 
tence of a mixed phase (MP) of quark and hadronic matter over a finite range 
inside the star. We have systematically studied the influence of the model 
parameters of the QP (bag constant B, coupling constant g) on the EOS (see 
Figs. |], ^), the phase transition densities (see Figs. ||, |7]), the mass-radius rela- 

4 For example the mass difference of the Hulse- Taylor pulsar PSR B1913+16 and its 
companion is AAf w 0.05 Mq while the difference of PSR B1534+12 and its companion is 



even of the order of only AM « 10 3 M© 14 1. Of course the two stars of a rising twin must 
not necessarily build a binary system. 
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tions of the stars (see Figs. ||, |) and the corresponding internal structure (see 
Figs. 11, [12]). This analysis provides us with the information which parame- 
ters lead to neutron stars with pure quark phase core (QC), mixed phase core 
(MC) or with no phase transition at all (HC). We found that - almost inde- 
pendent of the model parameters of the quark phase and the applied hadronic 
EOS - the radius of neutron stars possessing a MC can at most be about 20% 
more compact than a pure hadronic star of the same mass. A neutron star 
possessing a MC would therefore - at least from its radius - be hardly distin- 
guishable from a pure hadronic star. Stars with a QC, however, are typically 
20 — 30% more compact than a hadronic star. This is depicted in Fig. pi] and 
Figs. [H| |14[ The limitation to the upper value of 30% comes from our reason- 
able requirement that the deconfinement phase transition density should not 
appear below normal nuclear matter density. This approximately corresponds 
to i? 1 / 4 >165MeV (c.f. Figs |5], [7|). Within our model, the minimal radius of a 
neutron star is about 9 — 10 km at a corresponding mass of 1.4 — 1.5 M . Such 
small radii of 9 — 10 km are only reached in neutron stars with a pure QC (as 
large as 4 — 7 km) and could not be explained using the pure hadronic models 
applied here. On the other hand, our models fail to explain radii of R < 9 km. 
Although some radius estimates of neutron stars suggest quite small values 



H| [16], [17|, |18|, [Hj, ^T|, the experimental limits (e.g. of distance measurements) 
and the uncertainties in the interpretation of the available data (e.g. by black- 
body fits) do currently not allow for a definite confirmation. Nevertheless, a 
future experimental confirmation of extreme small radii would offer a unique 
possibility to place new stringent constraints on the EOS at high densities. 
It is therefore necessary to study the theoretical limits of the compactness of 
neutron stars with inclusion of different softening mechanisms to disentangle 
the various scenarios in view of more restrictive future radius estimates. 

In Sec. |5| we have discussed the theoretical possibility of a third family of 
compact stars which might exist besides the two known families of neutron 
stars and white dwarfs (c.f. Fig. ]T5j ). Within our model, we have shown that a 
deconfinement phase transition can explain the existence of such third family 
(c.f. Fig. [16]). Compared to stars of the neutron star family, the stars of the 
third family can have similar masses while their radii can be up to about 3 km 
smaller. Without refering to our particular model for the EOS we argue that 
the availability of mass and radius measurements of only two compact stars 
can reveal the existence of a third family if the larger star is the (slightly) 
more massive one. We refer to such pair of stars as a rising twin (c.f. Fig. [TJ). 
By means of Gerlachs criterion (Sec. |5.1|) the existence of a third family is 



equivalent to the existence of a phase transition in the EOS at large densities 
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(which not necessarily need to be the deconfinement phase transition). If 
reliable mass and radius measurements would be available, the detection of 
a rising twin could therefore serve as a novel signature for a phase transition 
inside neutron stars. 
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Appendix 

In this appendix we intend to show that a simplified treatment of the phase 
transition construction like the one which appears in a one-component system 
(instead of a two-component system which neutron star matter in fact is) 
leads to an overestimation of the phase transition density. This is worth to 
point out since within a simplified treatment which occassionaly is utilized for 
the construction of the phase transition one might exclude the possibility of 
quark phases inside neutron stars |fEZ| . In a correct two component treatment, 
however, the onset of a mixed phase (MP) or a quark phase (QP) might already 
be occured inside the star. Compared to the simplified treatment this results 
in more compact neutron star configurations due to a larger density range 
occupied by the soft MP or QP EOS. 

The two-component system of neutron star matter in weak equilibrium 
can be reduced to a one-component one by requiring charge neutrality of the 
hadronic phase and the quark phase independently (i.e. locally) and not glob- 
ally in the MP. As already discussed in Sec. |2.3| , this reduces the number of 
independent chemical potentials (components) from two (e.g. fj, e ) to one 
(e.g. /i„). (In the following one argument (p n ) denotes charge neutral prop- 
erties while two arguments (/i n ,/z e ) denotes charged ones.) The condition of 
local charge neutrality is of course unphysical in the sense that charge neutral- 
ity in the MP can also be achieved by means of the weaker (global) condition 
Eq. fljjp. Furthermore such a construction is thermodynamically incorrect if 
applied to neutron star matter. Whereas the Gibbs condition 

PHp(^n) = PQp(^n), (16) 

of a one-component system does ensure the mechanical equilibrium and the 
chemical equilibrium of neutrons (and therefore of all uncharged particles) one 
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finds that the chemical equilibrium of the second component (/x e ) is in general 
not fulfilled in the MP i.e. (cf. Fig. |) 

/if P (/in) ^/ifVn). (17) 

Consequently, all charged particles of neutron star matter do not fulfill the con- 
dition of chemical equilibrium. The Gibbs condition (|II]) leads to the familiar 
first-order phase transition with a constant pressure MP. Since we know from 
the equations of hydrostatic equilibrium - the Tolman-Oppenheimer-Volkoff 
equations [[01 - that the pressure has to increase if we go deeper into the star, 



a constant pressure MP is strictly excluded from the star. That is not the 
case for a two-component phase transition calculation, since the pressure is 
increasing with density even in the MP (see e.g. Fig. H). 

To show that the one-component treatment leads to an overestimation of 
the deconfinement phase transition density we go back to Fig. [| where /j$ 
marks the critical chemical potential necessary for a HP— >MP phase transition 
of the two-component system. All we have to show is that at //W the one- 
component system is still in its HP. In other words 

PHP(^ ] ) > PQP^ ] ) (18) 

must hold since at a fixed chemical potential the phase with the larger pressure 
is the physically realized one. The point CD in Fig. |2| and Fig. |3] is defined by 
the phase equilibrium between the charge neutral HP and the charged QP, 
which reads 

Php(^ ] ) = Pq p(^\^ ] ). (point [TJ in Fig. 2) (19) 

Now we know that at fixed neutron chemical potential the charge density of 
the QP is given (in units of |e|) by 

QP, \ 9p QP (fJ, n , fl e ) 

Pc(Pn,Pe)= ^ • (20) 

OPe 

Since the QP is negatively charged at point CD (and everywhere above the 
charge neutral QP curve in Fig. ^]) we find that the right hand side of (p0|) must 
also be negative. Therefore, at fixed fi n (say the pressure is decreasing if 
we decrease \i e from //W to the lower value at which the QP is charge neutral. 
This means 

Vqp{pI\^)>Vqp{pI\ (21) 
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Putting ([T^) and ([ZT]) together, we finally obtain (|T8|). This shows that due 
to the treatment of the phase transition as a simplified transition of a one- 
component system, the transition densities are shifted to higher values as 
compared to the ones expected from the correct treatment of neutron star 
matter as a two-component system. 
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